Gravitational Trapping of the Bulk Fields in 6D 
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We present two new solutions to Einstein's equations in (1 + 5)-spacetime with a positive bulk 
cosmological constant. One solution has increasing and another solution decreasing bounded scale 
function <f>(r) without singularities in the range from the origin r = to the radial infinity r = +oo. 
For the both solutions it is shown that all local fields are localized near the origin r=0 in the extra 
space through the gravitational interaction. 
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The idea that our observable 4-dimcnsional universe may be a brane extended in some higher-dimensional spacetimc 
with non-compact extra dimensions and non-factorizablc geometry has a long history [1] and has been the subject of 
many recent studies. Large extra dimensions offer an opportunity for a new solutions to old problems (smallness of 
cosmological constant, the origin of the hierarchy problem, the nature of flavor, etc.). In such theories it is assumed 
that all the matter fields are constrained to live on the 3-brane and the corrections to four dimensional Newton's 
gravity low from bulk gravitons are small for macroscopic scales. But this models still need some natural mechanism 
of localization of known particles on the brane. The brane solutions and questions of matter localization on the 
brancs in higher-dimensional bulk spacetimes have been investigated in various papers [2, 3, 4, 5]. In our opinion 
the localization mechanism must be universal for all types of 4-dimensional matter fields. In our world the gravity 
CT) i is known to be the unique interaction which has universal coupling with all matter fields. So, if extended extra 
s ! ' dimensions exist, it is natural to assume that trapping of matter on the brane has a gravitational nature. It is of great 
interest that recently in [G] it was found the model in the brane world where all the local bulk fields (ranging from the 
spin scalar field to the spin 2 gravitational field) are localized on the 3-brane only by the the gravity. The solution 
^ I , was found in (l+5)-spacetime for the positive bulk cosmological constant A > and has increasing scale factor 4>{r) 
ly-j ■ asymptotically approaching finite value at the radial infinity. In the paper [7] this solution was extended to the case 
, of a general (p-l)-brane model with codimension n in general D = p + n space-time dimensions. 

In our previous article [5] we have introduced the new 3-brane solution in (l+5)-spacetime which also localizes 
all the local bulk fields but in contrast to the solution found in [6] has decreasing scale factor 4>{r) asymptotically 
approaching finite non-zero value at the radial infinity. Such solution exists for the negative cosmological constant 
A < 0. In this article we present two new regular solutions (one with increasing and another with decreasing scale 
factor 4>(r)) to Einstein's equations in (1 + 5)-spacetime with a positive bulk cosmological constant A > 0, and for the 
both solutions we explicitly show that whole local fields (spin scalar field, spin 72 spinor field, spin 1 gauge field, 
spin I/2 gravitino field and spin 2 gravitational field), as well as totally antisymmetric tensor fields, are localized near 
(— 1 \ the origin r=0 in the extra space. In addition we investigate the technical reasons of the localization of all bulk fields. 

Let us begin with the details of our solution. In 6D the Einstein equations with a bulk cosmological constant A 
. £h ' an d stress-energy tensor Tab 

Rab - -g A BR = jjii^AB + T A b) (1) 
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can be derived from the whole action of the gravitating system 



4 



— {R + 2k)+L 



(2) 



Rab j R, M and L are respectively the Ricci tensor, the scalar curvature, the fundamental scale and the Lagrangian of 
matter fields (including brane). All of these physical quantities refer to (1 + 5)- space with signature (H — ■•• — ), capital 
Latin indices run over A,B,... = 0,1,2,3,5,6. Suppose that the equations (1) admit a solution that is consistent 
with four-dimensional Poincare invariance. Introducing for the extra dimensions the polar coordinates {r,0), where 
< r < +00 , < 9 < 2tt , the six-dimensional metric satisfying this ansatz we can choose in the form [3, 6] : 

ds 2 = <P 2 (r) nafi (/) dx a dx f3 - g (r) (dr 2 + r 2 d9 2 ) , (3) 

where small Greek indices a, (3, ... = 0, 1, 2, 3 numerate coordinates and physical quantities in four-dimensional space, 
the functions 4>(r) and g(r) depend only on r and are cylindrically symmetric in the extra-space, the metric signature 
is (H — ... — ). The function g(r) must be positive to fix the signature of the metric (3). 
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The source of the brane is described by a stress-energy tensor Tab also cylindrically symmetric in the extra-space. 
Its nonzero components we choose in the form 

T a p = -g a pF (r) , Tij = -g tj F (r) , (4) 

where we have introduced two source functions Fq and F, which depend only on the radial coordinate r. 

By using cylindrically symmetric metric ansatz (3) and stress-energy tensor (4), the Einstein equations become 

A" , A n , g" g' 2 , g' _ g ,„ ^ , g A phys 



,flV , A' _ 9 (F n J Khys 



6 ^ + 2 "gJ + ^ = 

A ' 2 J# 9 ,„ a \ , „ 9 K 



where the prime denotes differentiation dj dr. The constant A p h ys represents the physical four-dimensional cosmolog- 
ical constant, where 



B (4) n p(4) _ ^phys ( os 

H a0 ~ -j9otpn> - M 2 9af3- (») 
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In this equation R ^l, and Mp are four-dimensional physical quantities: Ricci tensor, scalar curvature and Planck 
scale. 

In the case h- p h ys = from the equations (5), (6) and (7) we can find 

F' + 4^ (F - F ) = 0, (9) 

<P 

g = , S = const, (10) 

r 

<t>" „ rg 

where 6 denotes the integration constant. The (9) represents the connection between source functions, it is simply 
a consequence of the conservation of the stress-energy tensor and can be also independently derived directly from 

d a t£ = 0. 

Let us take the cylindrically symmetric source functions in the form 

F (r) = h<t>- 1 ~ f2<t>-\ (12) 

foW^/if 1 -^- 2 , (13) 

where fx and are some positive constants. These source functions satisfy the equation (9). 
Taking the first integral of the last equation [3] , we get 

r4>= -WM-*V IA 0+ 3A,) V (14) 
where C is the integration constant. Then the equation (14) can be presented in the form 

r0 ' = ~iw ( ^ di)( ^ d2)+ 5' (i5) 
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where d\ and c?2 denote the roots of the quadratic equation cf> 2 — ^^<t> + §7^ = with respect to the variable 4>: 



bh 64A/ 2 \ 5/i / / 64A/ 2 , 



Let us suppose A > and 6 ^Y 2 < 1. Then for the di and <i 2 we have g? 2 > rfi > 
Setting C = 0, introducing the parameters 



0. 



a = Y^4' & = a(d 2 -di), (17) 
and imposing the boundary conditions for the solution in the form 

<f> ( r )\r=0 = C0n3t > °> <f> (0 U+oo = COnst > 0> ( 18 ) 



S(0lr=0 = COnsi ' ffWU+oo = COnst > ( 19 ) 



we can easily find two solutions of the equation (14) in the following cases: 
i) A > 0, b > 2, 5 > 



l + cr b ' JV ' _i/ (1 + cr 6 )* 



Hr)= \ ! b , 9(r) = Sb(d 2 -d 1 )- -5; (20) 



ii) A > 0, 6 < -2, (5 < 

^( r ) = 161 T . g(r) = \S\\b\(d 2 -d 1 )— (21) 

where c > is some positive integration constant. As we can see these solutions exists in the case of positive bulk 
cosmological constant A > 0. The scale factors </>(r) and g(r) of these solutions are smooth bounded functions of radial 
coordinate r. The scale factor <f>(r) of the first solution monotonically increases from the value equal to g?i > at the 
origin r = in the extra space and asymptotically approaches the value equal to d 2 > at the radial infinity r = +00, 
while the scale factor 4>(r) of the second solution monotonically decreases from the value equal to d 2 > at the origin 
r — in the extra space and asymptotically approaches the value equal to d\ > at the radial infinity r = +00. As 
we mention below one of the reasons of the localization of all local fields is connected with the properties of the scale 
factor g(r). Let us briefly consider the properties of this function for the first solution (20) (for the second solution 
(21) the results are the same) for various values of parameter b. In the case |6| = 2, at the origin r = and at the 
infinity r = +00 in the extra space this function has values equal to 2c<5(<i 2 — d±) and respectively. The derivative of 
this function g'(r) is negative in the range < r < +00 except at the origin of extra space r = where it is equal to 0. 
So in this case g(r) is monotonically decreasing function between r = and r = +00. In the case \b\ > 2, this function 
is equal to at the origin r = in the extra space, it increases in the range < r < r max , approaches the maximum 
value at the r max , then decreases in the range r max < r < +00 and monotonically approaches the asymptotic value 
equal to at the radial infinity. Here we have introduced r max which can be found from the equation g'(r) = and 
is given by 



6 + 2 



(22) 



In this case the derivative g'(r) is smooth bounded function except for the case \b\ < 3 when the derivative g'(r) 
approaches the asymptotic value equal to infinity only at the origin r = 0. 

From the formula (22)it is easy to notice, that in the case of fixed parameter b the r rnax tends to the origin r = 
in the extra space when the positive integration constant c tends to infinity. In the case c ^> 1 the scale factor g(r) 
has a (5-function-like behavior and therefore it can provide the sharp localization of bulk fields near the origin r = 
in the extra space (for the second solution (21) we have the same results in the case c«l). 

Now we turn our attention to the problem of the localization of the bulk fields on the brane in the background 
geometries (20) and (21). Of course, in due analysis, we will neglect the back-reaction on the geometry induced by 
the existence of the bulk fields, and from now on, without loss of generality, we shall take a flat metric on the brane. 
In the following we examine only the first solution (20), since the results for the second solution are the same (21). 
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We start with a massless, spin 0, real scalar coupled to gravity 

1 



So 

The corresponding equation of motion has the form 

1 



-.0 



M 



d 6 xv^gg AB d A <5>d B $. 



Zgg MN d N $) = 0. 



(23) 



(24) 



It turns out that in the background metric (20) the zero-mode solution of (24) is 3>o (x M ) = v (x^) po, where p = 
const, and v (x M ) satisfies the Klein-Gordon equation on the brane rj^d^d^v (x M ) = 0. Substituting this solution into 
the starting action (23), the action can be cast to 



So 



1 2 
nPo 



2tt 



dO 



dreffgr / d A xy/=rrif v d^v (x a ) d v v (x a ) + 



(25) 



(+oo) 



0(0) 



4> 2 d<t> / d 4 x^jn^ d^v (x a ) d v v {x a ) + ... 



(26) 



{4 - 4) J ^xyf^rpfdyf} (x a ) d v v (x a ) + 



(27) 



The integral over r in (25) is finite, so the 4-dimensional scalar field is localized at the origin r — in the extra space. 
For spin 1/2 fermion starting action is the Dirac action given by 



Si 



d b x^/~ Q g^iT M D M ^ 



from which the equation of motion is given by 

T M D M y> = (T^Djj + T r D r + T e D g ) * = . 



(28) 



(29) 



We introduce the vielbein h\ through the usual definition g^B — ^a^bVab wncre A,B, ... denote the local Lorentz 
indices. T A in a curved space-time is related to r y A by T A = /i~7 A . The spin connection co^ N in the covariant 
derivative Dm^ = [Om + \ u m N Imn) ^ 1S defined as 



MN 



\h NM (d M h% ~ d N hl) - \h™ (d M h% - d N h%) - \h PM h^ (d P h Q ~ R d Q h pR ) hi. 



After these conventions are set we can decompose the 6-dimensional spinor into the form <3> (x M ) = ij; (x M ) A (r) ^ e . 
We require that the four-dimensional part satisfies the massless equation of motion "f^d^ip (x' 3 ) = and the chiral 
condition j r i/j (x^) = ip (x^). As a result we obtain the following equation for the s-wave mode 



d r 



1 d r [rgi 

2 rg^ 



A{r) = 0. 



The solution to this equation reads: 



A (r) = Ao<p 2 g 



with Ao being an integration constant. Substituting this solution into the Dirac action (28) we have 

drift' 1 g* [ d^x^^4>i^d^ + ... = 



(30) 



(31) 



(32) 
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= 2vAl^5{d 2 -d 1 )c ^ ———— J d^x^P^d^ + ... . (33) 

In the case b > 2 the integral over r is obviously finite. Indeed, the integrand in (33) scales as r _ 3 fc_1 at the radial 
infinity and is the smooth functions between r = and r = +oo, so this integral over r is finite. Hence the spin I/2 
fermion is localized at the origin r = in the extra space only by the gravitational interaction. 
Now let us consider the action of U(l) vector field: 

Si=-\J d & x^~gg AB g MN F AM F BN , (34) 

where Fmn = 9m An — On Am as usual. From this action the equation of motion is given by 

-T=S=dM (^99 MN 9 RS F NS ) = 0. (35) 



By choosing the gauge condition Ag = and decomposing the vector field as 

Afj, (x M ) = 0|t (a") J2 °™ W eM > ( 36 ) 



A r (x M ) = a r °m (r) e iW , (37) 

l.m 

it is straightforward to see that there is the s-wave (I = 0) constant solution a m (f) = oq = const and a r = const. In 
deriving this solution we have used 9 M a M = <9 M / M „ = with the definition of / M „ = d^a v — d u a^. As in the previous 
cases, let us substitute this constant solution into the action (34). It turns out that the action is reduced to 

Si = -\ol drgr J d^Xy/^m^rjT Mp v + - = (38) 



/ dS / d^xV^n^t" UUu + ... = -^ral (d 2 - di) / d^xV^VV^Ufa + 
J4>{a) J * J 



(39) 



As we can see from (39) the integral over r in (38) is finite. Thus, the vector field is also localized. 

Next we consider spin ^/r, field (the gravitino). We begin with the action of the Rarita-Schwinger gravitino 



5*3 



(Fx* 



>g^ A iT [A T B Y c] D B ^ c , 



(40) 



from which the equation of motion is given by 



T [A T B T c] D B ^c = 



(41) 



Here the square bracket denotes the anti-symmetrization and the covariant derivative is defined with the affinc 



connection Tg C = /i~ ( d B h 



b"c 



, , BC h ~ 
uj b n cc 



by 



D A ^ B = d A ^B ~ ^ AB ^C + -uTl AB *B 



(42) 



After taking the gauge condition = wc look for the solutions of the form ^ (x A ) = {x v ) u (r) ^ e lW , 
\J> r (x A ) — ip r (x v ) u (r) J2 e%W where ^^(x") satisfies the following equations j u ipv = d M Vv = 7^ 7 P 7 r '^p'0r = 0, 
7 r Vv = ipu- For the s-wave solution and ip r {x u ) — we have for the equation of motion (41) the following form 



Or 



3 6' 



1 d r [rg-' 

2 rg5 



u(r) = 0. 



(43) 
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The solution to this equation is 

u (r) — u (j)~~i g~^r~? , (44) 
with Mo being an integration constant. Substituting this solution into the action (40) we get 

r+oo p 

S3=2ttuIJ dr<j)- 2 g^ J d i x^f^^i 1 [il Yl p] du^p + ■■■= (45) 

r+oo r^ b (l + cr b ) f 

= 2nu 2 ^/5(d 2 -d 1 )c / dr— -i > d *xV=Wp^l v l p] d u iP p + ... . (46) 

Jo r (di + d2cr b ) J 

As in the case for the spin I/2 fermion (33), when b > 2 the integrand in (46) scales as r~^ b_1 at the radial infinity 
and is smooth function between r = and r = +00, so the integral over r is finite. This means that spin 3 / 2 field is 
also localized at the origin r = in the extra space. 

Now let us consider spin 2 gravitational field. In this case we consider the spin-2 metric fluctuations H pv : 

ds 2 = {4> 2 (r) r] aP (x v ) + H aP } dx a dx< 3 - g (r) (dr 2 + r 2 d0 2 ) . (47) 
The corresponding equation of motion for the fluctuations has the following form: 

-4=9a (V^g AB d B H^) = 0. (48) 



/ e 9 

From the 4-dimensional point of view fluctuations are described by a tensor field H^ v which is transverse and traceless: 

0„ftf = O, H£ = 0. (49) 
This tensor is invariant under 4-dimensional general coordinate transformations. We look for solutions of the form 

H a0 = h a0 {x u ) r m (r) cxp(i6l), (50) 
m I 

where d 2 h fJiV (x a ) = m l h il y(x a ). It is easy to show that the equation of motion (48) has the zero- mass (mo = 0) and 
s-wave (I = 0) constant solution tq = const. Substitution of this zero mode into the Einstein-Hilbcrt action leads to 

r+OO r 

S 2 - 2ttt 2 J dr(/) 2 gr J d i x [d p h af3 d p h a , 3 + ...] = (51) 



r<t>(+°°) r 9 I" 

2ttSt 2 / <j> 2 d(j) / d 4 x [d p h a(3 d p h af s + ...] = -itSt^ (d\ - df) / d 4 x [d p h a/3 d p h af 3 + ...] . (52) 

J (l>(0) J 3 j 



The integral over r is finite, so the bulk graviton is localized at the origin r = in the extra space. With our 
background metric the general expression for the four-dimensional Planck scale Mp, expressed in terms of M, is 

r+ca r+00 n 

M 2 P = 2ttM 4 / dr(f> 2 gr = 2nSM 4 / dr(f> 2 (j)' = -TrSM 4 (dl - d\). (53) 
Jo Jo 3 

The inequality M <IC Mp is possible by adjusting M 4 and the product 8(d\ — df), and thus could lead to a solution 
of the gauge hierarchy problem. 

Finally, let us consider the totally antisymmetric tensor fields. The action of /c-rank totally antisymmetric tensor 
field Ak is of the form 

Sk = ~ j F k+1 A*F k+1 , (54) 

where F^+i = dA k . The corresponding equation of motion is given by 

dA*F fe+ i = 0. (55) 
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It is easy to show that A W(I2 i = a AtlA , 2 ... Mfc (x u )uq with u$ = const is a solution to the equation of motion (55) if 
d A */ = where / = da. Substituting this solution in the action (54) leads to the expression 

S k ~ / + dr^gr f f k+1 A */ fc+1 + ... = (56) 



= S y 4> 2 - 2k d<jy J h+i A */ fc+1 + ... = (3 _ 2fc) (^- 2fe - d\- 2k ) J f k+1 A + ... (57) 

As we can see from (57) the integral over r in (56) is finite, so the totally antisymmetric tensor fields are also localized 
by the gravitational interaction. 

In conclusion, in this article we have presented two new solutions to the Einstein's equations in the (1 + 5) -spacctimc 
with increasing and decreasing scale factors <f> (r) . These solutions arc found for the positive bulk cosmological constant 
A > 0. In addition for the our solutions we have presented a complete analysis of localization of a bulk fields at the 
origin in the extra space via only the gravitational interaction. There are two technical reasons of the localization 
of all bulk fields. The first reason lies in the fact that the scale factors (j>{r) are smooth bounded functions without 
singularities from the r = to the radial infinity. The second one, which in our opinion is the main reason, is connected 
with the function g(r). As it was mentioned above in the case c> 1 for the first solution (20) and in the case cC 1 
for the second solution (21)this function has a J-function-like behavior, and in contrast to the model considered in [ 
in our model the wave functions of the zero-mode solutions of the bulk fields along the extra dimensions are peaked 
at the location of the branc. 
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